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Abstract 



An extension is given for the standard two component model of 
adiabatic, Born-Oppenheimer (BO) electronic states in a polyatonic 
^ ■ molecule, by use of Mathieu functions of arbitrary order. The curl or 

>0 I compatibility conditions for the construction of a diabatic set of states 

based on a finite- dimensional subset of BO states are not satisfied 
exactly. It is shown, however, that, by successively adding higher or- 
(^ • der Mathieu functions to the BO set, the compatibility conditions are 

Tij" ' satisfied with increasingly better accuracy. We then generalize to sit- 

^^ . nations in which the nonadiabatic couplings (the dynamic corrections 

c/3 ' to the BO approximation) are small (though not necessarily zero) be- 

y i tween a finite-dimensional BO subset and the rest of the BO states. 

C/3 ' We prove that approximate diabatic sets exist, with an error that is of 

j^i the order of the square of the neglected nonadiabatic couplings. 

Ph' 

PACS number(s):31.15.-p, 31.30.-i, 31.30.Gs, 31.70.-f 

^^ . keywords: Born-Oppenheimer states, nonadiabatic coupling, Math- 
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1 Background and Introduction, preceded by a 
Homage 

Some of the authors in the article have had the good fortune that their 
pathways crossed with those of Per-Olov Lowdin, be this in Slater's Group 
at MIT, in Menton or at the Sanibel Workshops. They keep the memory of 
an endowed and innovative researcher, as well as of a Scientist Statesman, 
who followed in the footsteps of another great Scandinavian, Niels Bohr. 

The question of whether it is possible to effect in general an adiabatic 
to diabatic transformation (ADT) has been recently reopened [1],[2],[3]. 
One formulation of the issues [4] -[9] is whether a strictly diabatic basis in 
which linear combinations of a small number of Born-Oppenheimer elec- 
tronic states are chosen in such a way that the nonadiabatic coupling be- 
tween the transformed states vanishes exists. 

In view of possible ambiguities let us define the terms: 

Ck = Cki^e, Xr) (1) 

is a solution of the Born-Oppenheimer wave-equation for the electronic part 
(electronic coordinates Xg, to be suppressed at a later stage), involving (as 
parameters) the nuclear coordinate set {Xr}. For any given value of {Xr}, 
the Cfc form a complete set, known as the adiabatic set or Born-Oppenheimer 
electronic states, labeled by A; = l..cx). Let us consider only the finite set 
with k = 1, , n (where n is in practice 2 — 4). We shall call this the P-set and 
designate its complement the Q set. (This nomenclature follows previous 
practice to decompose the full Hilbert space into two disjoint subspaces, P 
and Q [10]. Then the projection operator I over the full Hilbert space can 
be written as the sum of two projection operators in the sub- spaces, namely 
I = P + Q ) The challenge has been to find a so-called diabatic set 

related to the P adiabatic-set alone, such that the nuclear derivative term 
in the nuclear part of the Born-Oppenheimer equation corresponding to S,m 
vanishes. Baer [11] showed that this can be achieved with an ADT matrix 
Aij{Xr) that satisfies 

for i,j, k in P and all nuclear coordinate indexes r. (A tensor-algebra nota- 
tion is used here, but supplementing the double-index summation convention 



with specification of the range of summands. The symbol r after a comma 
represents, as usual in tensor algebra, differentiation with respect to Xr). 
Then 

(k, m in P) and the nuclear Schrodinger equation contains no first-order 
derivatives of the nuclear co-factors of ^k ■ In equation (3) the following 
integrals over the electronic coordinates appear: 

4m =< CfelCm.r > (5) 

('•) ^ (r) (r.^ 

km mk v / 

The initial value conditions for A^m for some chosen initial {Xr = Xro} are 

Akm{{Xro}) = ^km (7) 

This is a convenient set of initial conditions, which can, however, be general- 
ized to cases when initially ^^m is non-zero for any k and m both inside P, 
or both inside Q [2]. Reference [11] also noted the following compatibility 
(curl) conditions as a requirement for the solution of equation (3) : 

Jr) _ is) ^ Jr) (s) _ (s) (r) (^. 

' km,s ' km,r ' kn ' nm ' kn ' nm V"7 

(k, m, n in P; all unequal r, s. A shorthand way of writing this equation 
in terms of many-dimensional vector matrices is curlr = —Tii.T.) These 
relations arise from differentiation of equation (3) with respect to Xg, a 
further differentiation with respect to X^ of another equation of the form 
equation (3) but involving s instead of r, and subtraction. The actual 
requirement of compatibility is that 

(Jr) _As) w (Jr-) (s) _ (s) (r)^^A (q.\ 

\'km,s 'km,r)^mh — {'kn ' nm ' kn 'nmJ^mh \i>J 

(k, m, n, h in P). Satisfaction of the relations in equation (8) is a sufficient 
condition for equation (9) to hold for all and any Amh and for the existence 
of solutions in equation (3) . 

Having now formulated the general issue, we work out in the next section 
a case, involving Mathieu functions, in which equation (8) is not satisfied. 
These functions have already been used in the context of the ADT matrix 
[12]. We show here how the compatibility condition can be better and better 
satisfied by successively enlarging in a systematic way the finite subset (the 



P subspace). Having illustrated with the Mathieu functions the case of 
small non-adiabatic coupling between P and Q subspaces, we then return 
in section 3 to the general theory and show what modifications are needed 
in the ADT matrix, generally, when the non-adiabatic coupling between P 
and Q are small, of order e. In section 5 we employ an analogy to show that, 
as in the case of the ADT matrix, frequently a large Hilbert-space is needed 
formally , but in practice a restricted space suffices. 

2 Mathieu Functions as Adiabatic States 

An electronic Schrodinger equation for the angular electronic coordinate 6 
and the polar nuclear coordinates q and (/) was written in Ref. [12], Eq. (1), 
as 

i-^^el-^ - G{q)cos{29 -4>)- u{q, 4>)]^{9\q, 0) = (10) 

where E^i is a characteristic electronic quantity, G{q) is a nuclear-electronic 
interaction coefficient, frequently assumed to be proportional to the nuclear 
coordinate q and u{q, (j)) is the eigen-energy of the solution (being part of 
the adiabatic potential for the nuclear motion). Equation (10) is recognized 
as the Mathieu differential equation, which has enjoyed a wide literature 
and a variety of notations [14]-[16]. We shall use the solutions given by 
these sources (thus differing from the methodology in Ref. [12], in which 
the solutions were derived), but shall retain the symbols for the parameters 
introduced in Ref. [12]. However, we shall make life easier for the reader 
by giving the relations between the parameters. In the relations, we indi- 
cate the literature sources by adding as subscripts the initials of the authors 
surnames. In the spirit of Born-Oppenheimer approximation the small pa- 
rameter in the theory is 
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^=Hww=-<1ml=-^F (11) 



and the adiabatic angular coordinate for the electronic motion relative to 
the nuclear one is 

z^e-^ (12) 

The last relation implies that all derivatives with respect to (j) can be 
replaced by -g times the derivative with respect to z. In the notation of 



Ref. [15], the two families of solutions of equation (10) that are of interest 
to us have the form: 

CO 

ce2n+i{z, -x) = J2 ^2r|l(-x)cos(2r + l)z (n = 0, 1, ..) (13) 

r=0 

oo 

se2n+iiz, -x) = J2 Bl';+l{-x)sini2r + l)z (n = 0, 1, ..) (14) 

r=0 

The functions are conventionally normalized to vr, so that the squares of 
the coefhcients sum up to 1. Thus, when we use equations (13 - 14) for 
the wave- functions ^ of equation (10) (that are normalized in the square 
to unity) each function has to be divided by ^/n. When we number the 
adiabatic electronic wave functions according to the order of their energy 
surfaces (as is commonly done in molecular physics) by {m'), we have the 
correspondences (for x > 0): 

^(3') ^ ce3{z,-x) ^(4,) ^ ^63(^,-3:) ^(5,) ^ ce5{z,-x) 
\/vr ' \/tt ' ^/7^ 



Solutions of equation (10) of the form ce2n, se2n have vanishing matrix 
elements with those in equation (15) and can be disregarded. For small 
|x|, the following are the leading terms in the expansions of the Mathieu 
functions for n = [15] : 



2 2 

X , , , X , X , , ^ X 



cei{z, -x) = (1 - Y^)cos{z) + g (1 - g)cos(3z) + Y^cos{5z) + ... (16) 



2 2 

X , , , X , X, , , X 



seiiz, —x) = (1 )sin(z) -\ — (1 H — )sin(3z) -\ sin(5z) + ... (17) 

Still for small |x|, the following are the leading terms in A2^^i{—x) and 
B^l'+li-x) [14, 15, 17]: for r > n: 



for r < n 



[(r -ny.{r + n + 1)!] 4 
{n — r)l . X 



n—r 



[(n-r)!(2n)!]^ 4' 



+ 0(x"-^) (19) 



Now the eigenvalues u of equation (10) are of the order of ^ " J — — 
[12, 15, 16] (this can be seen from the leading terms of equation (13) and 
equation (14) in which equation (18) and equation (19) were taken into 
consideration), therefore the n = pair is well separated energetically from 
the n > states that lie higher. Adopting the nomenclature of section 1, 
we shall call the two states in equation (13) and equation (14) with n = 
the P subset. For x = 0, P is doubly degenerate; for small |x| it is nearly 
degenerate. We look at the nonadiabatic coupling matrices r within this 
subset. Using the definitions shown in equation (5) for r = q and <j) (see 
also equation (15) in Ref. [12]) and the expansion in equation (16) and 
equation (17) we obtain: 

t[?1 = (20) 

and 

(4>) _ (1 + i) .on 

agreeing with equation(16) in Ref. [12]. Turning now to our equation (8) 
, we write the curl-term on the left hand side in the present curvilinear, 
plane-polar coordinate system as: 



dq dcj) ^ 32g2 2,2El 



22 



el 



where we have used equation (11) and equations (21,20) in the last result. 

However, evaluation of the right hand side of equation (8) , — r x r, 

inside the P-set (l',2') shows that this is zero. Clearly, the missing part 

(r) 

comes from the P-Q inter-set nonadiabatic coupling terms r}^. When we 
evaluate — r x r to the lowest order in x, we find that this exactly matches 
the curl- value and that the missing value come from the n = 1 functions ces 
and ses. This means that enlarging the P set from n = to n = 0, 1 ensures 
the compatibility conditions for the solution of the ADT matrix, correct to 
the order of E~^ . To satisfy the curl-conditions to higher powers in E^^ (or 
x) will require bringing in higher n values. It will be presently shown that 
each higher power of x^ requires one further n. 

Table 1 shows the leading values for small |x|, as obtained from equation 
(18) and equation (19) . 

We now turn to the nonadiabatic couplings T^'^\j,a) = —T^^'{a,J), 
where J is one of the n = and a is one of n > 0. Since we calculate these 
correct only to x, we find from table 1, upon recalling the power series expan- 
sions equation (16) and equation (17) for cei and sei and the orthogonality 
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Table 1: Some leading terms in the coefficients A and B for small \x\ 



of the trigonometric functions, that only n = 1 contributes nonvanishingly. 
We now draw up the list of the nonadiabatic couplings labeled, express them 
in terms of quantum mechanical bra-kets of the derivatives and give their 
values. We give those values of r in the expansions equation (13) and 
equation (14) , which contribute. 



^1'3' 
^1'4' 
"2'4' 
^2' 3' 



X 



X 

'Sq 



< cei Ibices >^ ■ 
-^ < ce2n+i\l£\sei >«^ 



(r 
(r 
(r 



0) 
0,1) 

1) 
0,11 



(23) 



From these we calculate the vectorial cross-product in the ground [P = 
(1', 2')] doublet 

{txt)i>2' = 5][r(«)(cei,a)rW(a,sei)-T('^)(cei,a)T(«)(a,sei)] (24) 

Collecting all matrix elements, after noting the antisymmetric character of 
the coupling matrices and the negative sign before second term in equation 
(24) , one finds finally 



(r X r) 



X 



1'2' 



32q^ 



(25) 



This agrees with the negative of the curl, as required. The contributing 
intermediate states are seen to be the lowest energy set from Q (n = 1), 
which is energetically immediately above the P (n = 0) set. This balances 
the x^ term in the curl, equation (22) . More generally, as one goes in cei, 
sei (or in <|>(i ),<|>(2 )) Iq higher approximations in x, curlTi'2' will contain 
higher powers of x^. (Recall that x oc -g— is the small parameter in the Born- 
Oppenheimer approximation.) We wish to show now that each consecutive 
term of order x^** in the expansion of cur/ri/2' will be balanced by a further 
set of the higher energy wave functions, precisely up to ce2s+i, se2s+i and 
none higher. (Since the method of computing the matrices is laborious, we 



do not calculate the higher order terms in the curl or in the vector product 
(r X r)i'2', only show from where the latter arise.) The proof relies on the 
substitution of the dominant term, given in equation (18) and equation 
(19) , into the expansions (13)-(14) for the Mathieu functions. We thus 
have approximately, for small |x|, 



oo 



ce2n+i(^,-x)«i^ai;fiixl"'"lcos(2r + l)z (n = 0,l,..) (26) 



r=0 



Se2n+l(^,-x)«^62n+l^|n-r|^.^2r+l)z (n = 0, 1, ..) (27) 

r=0 

where a and b are (in the leading order) numerical. Next we form the 
nonadiabatic coupling vector-matrix r, whose elements are proportional to: 



r) 1 °° 

< ceil — |ce2„+i >^ - E «2r+i«2r+iV " r|xl"--l+'- (28) 



-| o -| oo 

<Cei\- — \se2n+l >«-E«2r+1^2r|i(2r- + l)xl"-^l+'' (29) 

Q oz q ^_Q 

and two further matrix elements of similar form for coupling with sei. To 
compensate the term of order x'^^ in curlTi'2' by a similar term in (r x r)i'2', 
we require that 

\n — r\ + r = s (30) 

Two cases have to be considered: n > r and n < r. In the former, equation 
(30) gives n = s. In the latter case, one gets from equation (30) after 
transposing 

2r = n + s (31) 

If n were greater than s, then 

2r = n + s<2n (32) 

contradicting the assumption that n < r. Thus, in all cases, only Mathieu 
functions of up to order n = s need to be added to the intermediate set to 
compensate the x^* term in the curl; higher order functions start with higher 
powers of x. This leads to a hierarchical extension of the P set to ensure 
the compatibility condition to successively better accuracy in the ground, 
(l',2') set. 



3 The Projected ADT Matrix 

We return now to the general case embodied in equation (8) and equation 
(9). 

It has been estabhshed (Ref. [11], Appendix 1) that if in equation (8) 
the summation over n is extended to P + Q, i.e. the intermediate states 

(r) 

run over the complete set, then equation (8) is satisfied for r^^ defined by 
(5)-(6). 

We shall now consider cases in which all r^^ are small, (say) of order e, 
for either k or m being in the set Q complementary to P, though they can 
be large for both k and m being in P or in Q. Physical conditions for the 
existence of such situations have been noted before as either large energy 
gaps between P and Q ([5], section III. A, end), or a strong vibronic coupling 
localized in the neighborhood of some reference configuration ([5], section 
III.B), with which the Q-set has small overlap. In section 4 we consider 
yet another situation, where the smallness of e arises from the small ratio 
between electronic and nuclear masses. 

Our aim is to show that equation (3) and equation (4) continue to be 
approximately valid, even though equation (8) is not formally satisfied. Our 
demonstration should disarm objections to some practical uses of the ADT 
matrix. It also extends a previous proof [13], that the diabatic set exists 

(r) 

when the above considered r^^ are exactly zero, to those more frequently 
encountered situations, where the non adiabatic couplings are non-zero, but 
small. 

Notation: Henceforth we shall use Latin subscripts {e.g.,j,k,m, ..) for 
the total Hilbert space, capital Latin subscripts (e.g., J,K,M, .. ) for the P 
subset, and Greek subscripts (e.g.,a, /3, ..) for the Q subset. In this notation 
Tj^ are all of order e and small. We now write out equations (3) with initial 
conditions (7) for the full Hilbert space. Solutions of these exist, since the 
curl-condition in (8) is satisfied for the full Hilbert-space. Explicitly, (3) 
reads as follows: 

AjK,r + t%Amk + rjj^AaK = (33) 

for derivatives fully within P, and 

AaK,r + tJ^IAmK + T^a'pAfSK = (34) 

for the inter - PQ matrix element derivatives, and a further equation (which 
will not be of interest to us) for the intra-Q derivative matrix elements. 



Noting that at {Xr} = {Xro} the initial conditions (7) hold, we assume 
now that the matrix elements at a general point {Xr} close to it can be 
expanded in powers of {Xr — Xro} as a Taylor series. The condition for 

(r) 

this assumption is that the matrix elements r^^ have no singularities in the 
neighborhood considered [2]. The expansion takes the following form: 

oo -. M M 

^jk{{Xr}) = ^ l\n[X] ••• X! (^'"1 ~ ^ri0)---{Xr^ - Xr^o) 
N=0 ri=l rjv=l 

^jk,ri...rNi{-^ro})] (35) 

in which the coefficients contain the derivatives ^jfc,ri...rjv({^ro}) to the 
A^'th order {ri...rj\[ contains N similar or dissimilar symbols referring to the 
nuclear coordinates, and M is the number of nuclear coordinates). We shall 
find that all inter-PQ A^j are of order e, but the corrections to the intra-P 
AjK are merely of order e^. 

To start, we prove this for the first derivatives AjK,ri{Xro}) and 
AaK,r{{Xro}) respectively. These can be immediately evaluated from equa- 
tion (33) and equation (34) , with the following results: In equation (33) , at 
{X = XrO} the third term is zero by the initial conditions (7) Akmi{Xto}) = 
6km, and this gives the second term only, which has the same result as ig- 
noring the Q set, that is: 

AjKAiXto}) + T^ti{Xto})AMK(.{Xto}) = (36) 

AjKA{Xto})+r^fk{Xto}) = 0. 

In equation (34) the third term is similarly zero by the initial conditions, 
thus: 

AaKAiXto}) + T^jl{{Xto})AMK{{Xto}) = (37) 

AaKA{Xto})+r^jk{Xto}) = 

we conclude that AaK,r{{Xto}) = ~TaKi{-^to}) ^^ °f the same order of e, 
due to T^xii-^to}) being of this order. 

In the appendix we establish by mathematical induction that the cor- 
rections to all the inter-PQ set derivatives AaK,ri...rj^{{Xto}) are of order 
e, and that the corrections to the intra-P derivatives AjK,ri...rNi{-^to}) are 
of order e^. Subsequent derivatives are obtained recursively, rather than 
through integration, so that the validity of the procedure depends on the 
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convergence of the Taylor series (which is assumed to hold close to the initial 
point). 

Summarizing, subject to the assumption of small (though not neces- 
sarily zero) inter-PQ non-adiabatic coupling and analyticity of the matrix 
elements, the original procedure of Ref. [11] restricted to a finite dimensional 
set (P) is approximately valid, and correct to order e^, no matter that the 
compatibility conditions are not exactly satisfied. Earlier, it was shown, for 
the model discussed here, that the discrepancy in the compatibility condi- 
tions (8) is e^ (e.g., Ref. [5, 9]); our result establishes the magnitude of error 
in the solution. 

4 An Inductive Proof 

We now establish by mathematical induction that the corrections to all 
the inter-PQ set derivatives AaK,ri...rNi{-^to}) are of order e, and that the 
corrections to the intra-P derivatives ^j_K,ri...rjv({^to}) are of order e^. 

This result was already proven for N = 1. We assume that the result 
is correct for iV — 1 and show that the result is also correct for N. We can 
write: 

AaK,ri...rNi{^to}) = C^ri...r„_ir„+i...rjv^aii',r„ ({-^to}) (38) 

Using equation (34) this can be written also as: 

AaK,ri...rNi{Xto}) = -C^ri...r„_ir„+i...rjv('^QM ^MK + ^J/^ A/3K)i{Xto}) 

(39) 

Or also as: 

AaK,ri...rNi{Xto}) = -dri...r^_^r^^^...rN{'^aM A.MK){{Xto}) 

- dr,...r^_,,^^,..,,,{T'::;'^ApK){{Xto}) (40) 

The first term in the right hand side of the above equation is proportional to 
''"aM ^^*^ ^^^ derivatives, and thus by assumption is of order e. The second 
term in the right hand side of the above equation contains derivatives up 
to order A^ — 1 of Ap^ and thus by the assumption of the induction is also 
of order e. Thus we established that the right hand side of equation (40) 
is of order e and also ^Q_ft-,ri...rjv({^«o}) is of order e. We conclude that 
AaK,rx...rN{{-XtQ}) IS of Order e for all A'' by induction. 

Next we want to establish that correction to ^j_ft:,ri...rjv({^to}) are of 
order e^. We established this result for N = \ inside the text. We proceed 
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by assuming that this result is true for A^ — 1 and wih prove its correctness 
for N, again we write: 

AjK,ri...rNi{^to}) =dr^,„r^_^r^^^...rNAjK,r^{{Xto}) (41) 

Using equation (33) we obtain: 

AjK,ri...rNi{Xto}) = -9ri...r„_ir™+i...rjv(Tj^ AmK + Tj^ A/SK)i{Xto}) 

(42) 
Or also as: 

AjK,ri...rNi{Xto}) = -'9ri...r„_ir„+i...rjv(rjj{^ AMK)i{Xto}) 

- dr,...r^_,r^^,...rAr%"'^ApK)i{Xto}) (43) 

The first term in the right hand side of the above equation contains intra-P 
terms which are not considered corrections. The second term in the right 
hand side of the above equation, which is considered a correction term, 
contains derivatives up to order A^ — 1 of A^k , which are of order e according 

(r ) 

to the theorem established above. When multiplied by derivatives of t^^ 
which are of order e by our assumptions, the correction terms obtained is 
thus of order e^ for N. Thus the result holds for any N, which is what we 
have set out to prove. 

5 An Analogy 

Very frequently one has a complicated Hamiltonian H, which has non- zero 
matrix elements Hnm between the components (designated n and m) in an 
extended Hilbert space. For many practical purposes, e.g. for numerical 
work, one uses a restricted, finite set (belonging to the P subset of the 
Hilbert space), such that P-Q matrix elements of the Hamiltonian are small. 
The procedure of solving the Schrodinger equation in this manner can be 
formally written as: 

HP = EpP (44) 

where, as before, P is a projection operator and Ep represents the set 
of eigenvalues corresponding to the set of solutions thus obtained. Pre- 
multiplication with Q, gives 

QHP = EpQP = (45) 
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since P and Q are disjoint. This implies that the procedure requires, for 
consistency, that QHP = 0, or that there be no matrix elements Hnm of 
the Hamiltonian connecting the Q and P sets. While this is, in general, a 
true requirement for exact solutions, in practice, when one seeks meaning- 
ful approximate results, one is content with the smallness of these matrix 
elements. We have here a further example of a rigorous condition being 
inconsistent with an approximate method. 

6 Conclusion 

In section 2 we have seen that for adiabatic electronic wave functions in 
the form of Mathieu functions the compatibility (or curl) conditions are not 
satisfied but, by successively enlarging the finite set, the conditions can be 
satisfied to increasingly better accuracy. As previously shown in (e.g.) [5] 
Eq. (2.54), a non-zero curl gives rise to a "gauge field tensor, the non- 
Abelian analog of the electromagnetic field" , also called the Yang-Mills field 
[18]. The use of Mathieu functions will thus enable a systematic study of 
the properties of this field in a molecular physical context [19]. 

We then discussed the question of the existence of a finite sized ADT 
matrix A (which is distinct from the question of its uniqueness or single val- 
uedness), being the solution of a finite number of partial differential equa- 
tions that fail to satisfy compatibility conditions. We have started with a 
complete set, for which the existence of solution is assured. Under condi- 
tions that the nonadiabatic coupling with the exterior of the set is small, we 
have truncated the size of the set and removed part of the finite sized ADT 
matrix A. The compatibility conditions are not met within the finite set, 
but the remainder of A satisfies the original differential equation approxi- 
mately, namely to an accuracy which goes as the square of the neglected 
nonadiabatic coupling. We conclude that the failure of the compatibility 
conditions yet permits the existence of an approximate ADT matrix (sub- 
ject to smallness of some nonadiabatic coupling terms). 
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